ENERGY DETECTION TROPOSCATTER SYSTEM WITH 

INBAND DIVERSITY 


By 

Gapt. M. JAYAKRISHNAN 




ENERGY DETECTION TROPOSCATTER SYSTEM WITH 

INBAND DIVERSITY 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 
MASTER OF TECHNOLOGY 



By 

Capt. M. JAYAKRISHNAK 


eBMMKTiiatT OF ejBgimaMb 

Diimti of rionimoGY, kanpuk 



TJ T •• '^'PUR 


te> Nu. i . 

IB OCT *-78 


551 ^ 


(^-D^^y-£N6 



ii 


GEETIFICATE 

This is to certify that the thesis entitled 'Ener^ 
Detection Troposcatter System with Inland Diversity* is a 
record of the work carried out by Capt. M. Jayakrishnan under 
my supervision and that it has not been submitted elsewhere for 
a degree. 


X*?) Septojiber, 1978 




Dr. P.K. Chatterji 
Assistant Professor 
Department of Electrical Engineering 
Indian Institute of Technology 
Kanpur. 


" 1 ^ . K A >» \ I t ( b !< 

* tip'fo'vf'd 

‘ di ({ » I tl if ! rt r i -/ 

^ ^ '■'-*< ( A'J i , i. I 

f ^ n t M ■ i d i K r v‘ 1 1 h f I'if* 

nj tae I/, fian jt 

Institute of 1 echu4)iugy K^/ipiK/! 
i).i(ed, l^L 



IIX 


ACKITOWLED GEMEIITS 

I am greatly indebted to Dr. P.Ii. CSaatterji for his 
constant and valuable guidance in the preparation of this thesis. 
He has evinced keen interest at every stage and helped me to 
complete this thesis. 

I would like to gratefully acknowledge all the help given 
by Dr. 7.P. Sinha and Dr. M. S. Zrishnamtirthy at various stages 
during the preparation of this thesis . 

I am also deeply obligated to Dr. V. Pajaraman, 

Dr. R. Sankar and Dr. P.R.Z. Sao who taught me during my stay 
here. 

I am singularly indebted to Lt.Col. Gurudev Singh for his 
valuable assistance while preparing and running the many computer 
programmes involved in this thesis, I have also drawn freely upon 
the skills of my colleagues and friends in preparing this thesis. 
Since the list of these becomes too extensive to permit individiaal 
acknowledgements, let each of them accept my sincere thanks for 
helping me In this effort. 

Fi na l l y, I also acknowledge with pleasure the excellent 
typing by Ifr. C.M. Abraham, 



LIST OF COilTENTS 


iT 

Chapter Page 

1. IHTRODUCTIOif 1 

1.1 Peculiarities of a Mobile Troposcatter 

System 1 

1.2 Diversity Techniques Suitable for Mobile 5 

Tropo 

1.5 Energy Detection 5 

1.4 Scope of the Work 7 

2. GHANKEL CHARACTERIZATIOE 8 

2.1 Introduction 8 

2.2 Gaussian 7SSUS Model 9 

2.5 Delay Power Spectrum 10 

2.4 The Tapped Delay Line Channel Model 11 

2.5 Calculation of the Delay Power Spectrum 15 

2.6 Parameters Used for the Simulation 15 

5. DIGITAL FILTERING AND PSEUDO RANDOM NOISE 19 

GENERATION 

5.1 Introduction 19 

5.2 Design of Digital Filters 19 

5.3 Butterworth Filters 20 

5.4 Design of Complex Digital Filter 21 

5.5 Details of Filters Designed 24 

5.6' Pseudo Random Noise Generator 24 

5.7 Generation of Uniform Distribution 25 

5.8 Generation of Gaussian Distribution 26 

4. SYSTEM DESCRIPTION AND DIVERSITY 27 

4.1 Signal Flow Model 27 

4.2 Simulated System 27 

4.3 Channel Foimiulation 27 

4.4 Modulation ^ 

4-5 Digital Filter 30 



Ghapt er 


V 


4.6 Gaussian Hoise Generation 30 

4.7 Energy Calculation 30 

4.8 Decision 30 

4.9 Computation of Statistics 31 

-,.10 Diversity 31 

4.11 Diversity Combining 31 

4.12 Signalling Scheme 32 

4.13 Signal Format Used 33 

5. SIMUMION AND RESULTS 39 

5.1 Introduction 39 

5.2 Main Program 39 

5.3 Results 43 

6 . CONCLUSION 47 

REFERENCES 51 



ABSTRAOT 


Reliable communications over fading dispersive channels 
requires the utilization of some form of diversity. Possible 
diversity techniques suitable for mobile troposcatter systems 
are examined, with particular emphasis on inband diversity. 
Some inband frequency diversity techniques are simulated and 
error probabilities are obtained for various SHRs. These 
error probabilities are usually less than those obtained by 
Bello [3]. 



CHAPTER 1 


IHTRODUCTIOH 

In racent years, a tactical troposcatter system on 
wheels, is a very attractive proposition for reliable communi- 
cations. This is especially true, when satellite communication 
for military usage on a large scale becomes prohibitive in 
cost. A mobile troposcatter ^stem is characterized by the 
constraints that it imposes on the size of the system. A tro- 
poscatter channel is characterized by the presence of fading 
and multipath besides the usual additive noise. Hence relia- 
ble communication over such fading dispersive channels requires 
the utilization of some form of diversity. 

1.1 Peculiarities of a Mobile Troposcatter System 

A mobile troposcatter system imposes certain constraints 
due to the’ mobile role it has to play. Mobility in this case 
implies ease of transportability, A mobile troposcatter 
system should be capable of being transported conveniently 
and the establishment of a communication link should be 
possible within a short time. These two basic requirements 
of a mobile troposcatter system effects many aspects in the 
design of the system. The first, involves a reduction in 
the size and bulk of the equipment. Therefore, the design 
should bepreferably be based on a single transmitter and 
receiver usirg only one antenna. Since the power plant has 
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also got to be mobile, there is a limitation on the amount of 
power available. The utilization of conventional forms of 
diversity like space and frequency involves increasing the 
size of tLi) system to accommocate the dual aitennas as in 
space diversity or the dual transmitters as in frequency 
diversity. These are, therefore, not advantageous in a 
mobile troposcatter system. The second requirement of esta- 
blishment of the link in minimum possible time results in 
limited alignment time. This is particularly so in conventiona.. 
space and frequency diversity techniques where there are more 
th^ one antenna to be aligned, 

N. 

A consistently higher degree of message recovery always 
involves some form of diversity technique , Hence, in the 
case of a mobile troposcatter system we therefore need some 
technique of diversity that does not have any of the constraints 
mentioned earlier. Some of the possibilities are ; 

a) Angle diversity 

b) Pclarization diversity 

c) Time diversity 

d) Multipath diversity 

e) In band diversity including modified forms of 
frequency diversity. 

These are d^J.t with in greater detail in the next section. 
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1.2 Diversity techniques Suitable for Mobile Tropo 
1.2.1 Angle Diversity 

Angle diversity can "be considered as a form of space 
diversity. While space diversity is obtained by the relative 
phasing among rays from one portion of the scatter volume, 
in angle diversity a number of transmitting beams are used 
each oriented towards a non overlapping portion of the scatter 
voliune and a similar number of receiving beams associated with 
each transmitting beam, This does not require any additional 
antennas, though a n-umber of feed cones are required at each 
antenna. Additional transmitter power is necessary. Angle 
diversity is considered as one of the best ways of achieving 
diversity nowadays in the context of mobile tropo systems [9], 

1-2.2 Polarization Diversity 

In this case we transmit signals polarized in one plane. 
In the channel depolarization occurs. Ihe signal is received 
independently on two orthogonal polarizations. The two 
received signals do not fade in a correlated manner. Since 
there are only two orthogonal polarization modes, this is 
essentially a dual diversity tedbnique. However, this 
technique is not useful in troposcatter since there is hardly 
any change in the polarization while the wave proj^gates. 
However, since the polarization does not change while propa- 
gating it i# pepsihl© to consider tran^dttipg, oi^ogophl, 
and receiving ,th« separately' |43» ■ ■ 
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1.2.3 Time Diversity 

In time diversity, the same information is repetitively 
transmitted at time intervals well separated as compared to 
the recipx’ocal of the average fading rate. Since the channel 
is time selective, the fade levels associated with the various 
repetitions will essentially be independent. Hence, appropriate 
combinations of the repetition will give a diversity performance. 
Unlike techniques like space diversity, time diversity is 
achieved without any increase to the total average received 
power per bit. 

1.2.4 Multipath Diversity 

This is another form of time diversity that can be achi- 
eved in situations where sufficiently wideband signals are 
transmitted and the multipath may be viewed as creating a 
series of echoes. The contributions ariving via the different 
parts of the multipath structure can be separated in the course 
of reception and combined in divei-sity. The Rake system uses 
such a technique. Here the . eceived signal plus noise is 
cross correlated against a replica of each possible transmi- 
tted waveform, at every possible delay and the results 
sumB^d with a maximal ratio combiner weighting using the 
estimate of the channel at that relative delay [4]- 

1.2.5 Ihband Diversity 

iB the diversity implicit in the channel du© to the 
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time and frequency selective fading that the signal undergoes 
in the channel. This is possible due to' the fact that 
pulses at different frequencies fluctuate in an uncorrelated 
manner. Thus if we transmit the same information at different 
frequencies, we have an inband frequency diversity. '[Jnlike the 
conventional frequency diversity technique where the seperation 
in frequency is larger, here we use a small seperation in 
frequency to achieve the effect of diversity. 

1.3 Energy Detection 

In order to obtain high data rates and reliable transmi- 
ssion on fading dispersive channels like troposcatter links, 
clever use must be made of the diversity inherent in the 
available transmission bandwidth. The basic approach is to 
transmit pulses whose time duration and/or bandwidth cover 
many selective fades in time and/or frequency and to use 
energy detection for pulse discrimination [3], 

It is possible to have both time and inband frequency 
diversity using energy detection. However, since the pulse 
width is very small compared to the fading time constant, 
time diversity is unsuitable as a large duration has to 
elapse to get two uncorrelated pulses. This means that the 
timie between two uncorrelated pulses is very large and this 
results in very low data rates. 
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However, if we use different frequencies in sniall time 
slots it is possible to acihieve frequency diversity without 
much loss in data rate. Here the bandwidth of the pulses 
is larger than the coherent bandwidth of the channel. The use 
of time-frequency pattern is one such technique. Here the 
spectral energy in the vicinity of the mark frequency is 
identified as mark and the energy around the space frequency 
as Space. 

Consider a non coherent PSE case with multiplicative 

fading occxiring at a rate faster or comparable to the pulse 

IKl. »vffvvU‘»W inAoflk IQ. 

rate. If T is the symbol duration and T^^jthe auto covariance 
of the complex Gaussian fading process, then over each S3n3ibol 
d\iration we have 

M = (1,1) 

a 

independently fading states occuring-. Bach of these M succe- 
ssive outputs spaced at intervals of T^ are independent complex 
zero mean Gaussian variates, each having a signal and 
a noise component. How if we have filters centered around mark 
and space frequencies and the samples of each filter square 
law combined over the symbol duration T, we can arrive at a 
decision based on the larger of the combined outputs. This 
in brief is the principle of energy detection. 

l^Ie it is passible to use other forms of diversity like 
angle diversity and polggrization diversity in a mobile 
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troposcatter system, one of the attractive methods is utili- 
zing inbond diversity using energy detection. This does 
entail a loss in data rate due to repeating the same informa- 
tion within a fixed bandwidth allocation. However, it should 
be possible to achieve hi^ data rates using good coding and 
decoding techniques, 

1*4 Scope of the Work 

The simulations reported here pertain to the study of 
inband frequency diversity using energy detection. We use a 
channel model as described in (2] for this simulation. The 
parameters of the link are as would be foxind in the case of a 
mobile troposcatter link and is described later. The channel 
characterization is discussed in Chapter 2, In Chapter 3, we 
discuss the complex digital filtering techniques used here 
aj3d the generation of complex Gaussian samples. In Chapter 4, 
the energy detection system and the different signalling formats 
for inband diversity are discussed. We have simulated three 
different time-frequency formats and the relative advantages of 
each are dealt with here. Chapter 5 deals with the simulation. 
The flow-chart is discussed here in detail, and results 
are disclosed. In Chapter 6, we conclude this thesis. 
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CHAPTER 2 

CHARBEL GSARACTERIZATIOB 


2,1 Introduction 

All radio tr.an.smission media may be considered as a 
randomly time variant linear channel. In a radio channel, 
like in any linear system, the transfer function or the impulse 
response provides necessary information for characterizing its 
input-output behaviour. Since the system functions are random 
processes that are at most quasi-stat ionary , such a characteri- 
zation is complex. One such idealized statistical model for 
representing fading dispersive radio channel is the (kiussian 
wide sense stationary uncorrelated scattering (WSSl®) model [1]. 

Certain simplifications can be made to the general chara- 
cterization of the randomly time variant channel, while trans- 
mitting digital signals. The simplifications become possible 
when the time and frequency selective behaviour of the channel 
may be regarded as wide sense stationary for time intervals much 
larger than the durations of the signalling waveforms that are 
being used. Such a situation arises when the channel contains 
very slow fluctuations superimposed upon more rapid flu- 
ctuations, which exhibit stationarity properties. 
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2. 2 Gaussian WSSUS Model 

The defining characteristics of a Gaussian wide sense 
stationary uncorrelated scattering model may be expressed as 
below ; 

The time variant transfer f^mction is a wide sense 
stationary complex Gaussian process in both time and frequency 
variables of the transfer function. That is, if we denote the 
time variant transfer function by T(f,t), the WSSUS property of 
the definition becomes 

E[T*(f,t) T{f+^, t+T)] = {2.1) 

where is the time frequency correlation fvinction. 

The time variant transfer function T(f,t) is the 
Fourier transform of the time varying imp-ulse response 
g(t, %). 

T(f,t) = (2.2) 

g(t,p = J df (2.5) 

The magnitude of T(f,t) has a Raleigh distribution function 
and the phase of T(f,t) is uniformly distributed [2], 

From the time frequency correlation function R(nu ,- 5 ) we 

get , 

a) the frequency correlation ftnaction 

(2.4) 
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b) the time correlation functions 
p(x) = R(0 ,t) 


(2,5) 


Since the channel has uncor related scattering from different 
paths, we have 


E[g*(t,’f) g(t+T,Ti)J = Qix,'l) 5(n-^) 


( 2 . 6 ) 


where 6(*) is the imit impxilse response function and Q(t:,^) 
is proportional to the o,uto correlation of the gain functions 
for paths providing delays in the range (^. , ? + in a 
tapped delay line model. This function is called the path 
gain correlation function or tap gain correlation function 
[ 2 ]. 


0(1, t) and R(-n., t) are Fourier transform pairs, 

(2.7) 

Q(t,^) = jR(n.,'T;)e^^’^^^-d-n- (2.8) 

2*3 Delay Power Spectrum 

The delay power spectrum QC^) as mentioned earlier is 
the Fourier transform of the frequency correlation function 

Kfn-). 

Q(t) = Q{0,53 = (2.9) 

it Ss etml to the power oontainet in the cfeamel tmasfer 
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function. In other words, the delay power spectrum Q('^) is 
the intensity of scattering for the delay . 

Physically it is the cross correlation between two 
carriers that are seperated by-OaEz and is a measure' of the 
channel^bandwidth. This correlation bandwidth becomes mono- 
tonically narrower for a given antenna size as the antenna 
separation increases, while for a given antenna seperation the 
correlation bandwidth increases as the size of the antenna 
increases [4]. From the width of Q(^) we can determine the 
amount of intersymbol interference and the extent of frequency 
selectivity. This width is the multipa.th spread of the channel, 

2.4 The Tapped Delay line Channel Model 

The ¥SSUS channel can be modelled as a tapped delay line 
with taps taken at intervals of ^ secs. [1], [4]> where 
is the sampling frequency. 

Assuming that the channel responds to frequencies in the 
interval 

< f <; fi+ vjz, 

then, we have 

-o 

g(t,p = 211 g(t,n/¥^) •ezp[a27if^(^ - n/¥^)] x 

stoo[ { 1 - 


( 2 , 10 ) 
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given the sample values of g(t,J ) at time intervals l/W^ 
apart, where 


Sine y = ^ 


If z(t) a nd w(t) are the complex channel input and 
output respectively, than we have by convolution. 




w(t) = / z(t-p g(t,v) d'e 
This results in 


( 2 . 11 ) 


w(t) =fiJl g(t, §7) . exp[3 2itfj_(;^- §7] z 

X Sinc[Wj_(’^~ I7]. z(t-J)] 

«o 

" - It)] . Sinci;¥^(5- | 7 ) 3 ]dJ 

"*" (2a2) 


i. 3Q. 

where gj^(t) = ^ g(t, 

The convolution of the input with the time function 
exp [3 2Ti:f j|^(t-:i/¥^)] Sine ¥^(t--n/V^)] 

corresponds to passing the input signal throiigh a band- 
limited filter [4], 

rect[(f-fj^)/W^] ezp|;-a2i£nf/¥^) 

For iiie oafee of WSSir0 channel with input frequency 
Constraints, the; input band liiaiting filter can be incorpcacated 
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into the taps, with the result that the correlation properties 
can he shown to he [3] 




= Sin<:[Wj^(J- |-)] -SinoiWj^C^- |-) ]Q(t, >)dj 




(2.13) 


where Q(t,^) is as defined earlier. 

If Q(t,^ ) changes little for changes in ^ 
of l/W, we have the approximation 



in the order 


m n 

(2.14) 

n = n 


Thus the tap ga.ins are uncorrelated random variables whose 
variance are proportional to Q(n/’'f^), the delay power density 
spectrum at the time n/¥^. Figure 2,1 shows the tax^ped delay 
line model. 


2.5 Calculation of the Delay Power Spectrum 121 

Figure 2,2 gives the geometry for calculiation of delay 
power spectmim. The details are as under : 

r = height of transmitting anten n a above the sea level 
s = height of receiving antenna above sea level 
d * straight line distance between transmitter and 


receive 
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R = radius of earth to sea level 
o 

= angle the radio horizon at the tr-iismitter makes 
with the line connecting the transmitter and the 
receiver 

= angle the radio hoi’izon at the receiver makes with 
the line connecting the transmitter and the 
receiver 

R = distance from the transmitter to differential 
scattering volume 

S = distance from the receiver to the dif f oi^ential 
scattering volume dV 
y = angle "between R and D 
0 = angle between S and D 

Q = scattering angle i.e. angle throiigh which 

incoming wave must be tiirned to sca.tter towards the 
receiver 

h = height of volume dV above line connecting trans- 
mitter and receiver 
c = velocity of light 

Using the bistatic radar equation 

dp/v (2.15) 

R2S®hO„ 

m 

Bello has derived in [2] that 
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V^/0 


Q(6)ci 


l+m7^ 


° <^{x^26 - '^ ) H(V26/x - 0^) 

o --2_ az (2.16) 


(x + X 


fjm 


whore the nornalized po,tli delay 




= path delay corresponding: to different i-o.! scattering 


R n 

volume = — - — , and 


X = a norma 


alized variable = V^/V^, 


m = scattering angle exponent. 


Por a Gaussian 




G(^) = H(<^) = exp(-^2/ar^) 


(2.17) 


where o may be expressed in terns of the half-powor bandwidth 
of the antenna as 


a = O.Sip ^ 

and = D/2R^ 

2.6 Parameters Used for the Simulation 


(2.18) 

(2.19) 


The simulation was conducted using the following 
parameters ; 

Distance between the transmitter and receiver = 175 miles 
ilatenna diameter = 28 feet 
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Half power beam width of the antenna 


IP 3andwif’xji 




HP Bandwidth 

Delay power spectrum time 


3 degrees 
70 mz 
12 Mz 
.6 fisec 


Humber of ma.rk space pairs 

Type - I 
Type - II 
Type - III 

Pulse width 

Type - I 
Type - II 
Type - III 

Humber of taps on tapped delay line 


= 1 
= 2 
= 1 

= .667 usee 
= 1 jisec 
= 1 |isec 
= 25 


A sampling freq,uency of 40 Miz was used for generating 
the tapped delay line channel. 



17 
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Fig, 2,1 . The Tapped Delay Line Model 
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CHAPTER 3 

DIGITAL EILTERIIIG MD PSEUDO RAHDOM ROISE 
GEHERATIOIT 


3.1 Introduction 

Digital filtering is the process of spectrun shaping. 

Even though the purpose of digital filtering is the same as 
that of analog filtering, tho realization of a digital filter 
is quite different from that of an analog filter. Theory of 
analog filters is based on mathematics of linear differential 
equations vrhile digital filtering is based on linear difference 
equations . 

3.2 Design of Digital Filters 

Digital filters can be designed by various techniques 
some of which are 

a) impulse invariant method 

b) bilinear transformation 

3.2,1 The Impuls e Invarian ce Method 

The impulse response k(t) of a continuous filter is 
defined as the inverse Laplace transform of its system fun- 
ction Y(s). 

k(t) = 



20 


Similarly tiie impulse response li(nt) of a filter 

is defined as the inverse Z tr-'.nsform of its sjrstem function 
H(z). 

Thus the condition that the impulse response of the 
digital filter is equal to the sampled impulse response of a 
cont inuous f ilt er 


l(nt) = k(t) t = ...... 

1 leads us to a digital filter. 

3.2.2 Bilinear Transformation Metho d 

If we have a stable analog filter defined by H{s), the 
frequency response of this filter can be found by evaluating 
H(s) at points on the imaginary axis of the s plane. Now if 
we ha,ve s replaced by a function of Z such that the imaginary 
axis of s maps onto the unit circle, then the resulting H(z) 
evaluated on the unit circle will take the same values that 
H(s) takes on the imaginary axis. The simplest rational fun- 
ction that maps the jw axis onto the unit circle is 

3.3 Butterworth kilters 

The class of analog filters with the system function 

= i -— — 

TT ) 


{3*2) 
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h-ving m poles at finite values of s and .a order zero for 
infinite s are called Butterworth filters. 

A third order Butterwortia Low pass analog filter has a 
system function [6J 

~ "(s+Sj^) (S+S 2 ) (s+s^) (3.3) 

with , 

Sg = 4-(l+3f3)w^, 

S 3 = i(l-3f3)w^ 

where w is the cut-off frequency, defined by YCjw ) = 0.707. 

c c 

The Z transform of such a filter is 


H(z) 


0Z~^ + DZ~^ 

~ 1 -w'„t7"^ 

(1-e ® Z^-^-Xl-Ze ° (Cos 



w^T)Z~^ + Z”^) 


(3.4) 


with 

-W T -W^T/2 t ■if-z 

C = w^ [e ^ + e Sin-^ w^T - Cos ^^5!)] 

-w T -3w^T/2 . 

D = w^[e ° - e ° Sin ^ w^T + Cos ^ wJT)] 


3.4 Design of Complex Digital filter [Yl 

Consider a bandpass function centered around w^. The 
real part r(t) can be defini^ as 
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r(t) = Re [z(t)] 

t 

= Re [e(t)e ] 


= Re [e(t) Cos w^t - Iia^eCt)] Sin w^t (3.5) 

This is a compact representation of the in phase and g.uadratare 

components of the bandpass fxmction. The effect of the modula- 
jw t 

tion by e is merely a translo,tion of the Bpectr\nn by an 

amount w . 

c 

Hence, in developing a complex filter we start with a reo-l 
filter centred n.round w = 0 and shift the centre point so that 
the desired amplitude and phase are now obtained rela,tive to a 
new frequency w., without introducing a conjuga.te filter at 
-w^, as would have occured in the case of a lojpass to band- 
pass transformation in a, real filter. 

Mathematically this shift is accomplished by replacing s 

by s - jw in the real filter transfer function H(s). 'Uhis 
° 3W t 

is equivalent to multiplying the impulse response by e 

The mapping 

Z = e®"^ (3.6) 


accomplishes this s plane transfer along jw axis, to a z plane 
shift for the digital filter [7]. 


Hence, we get 

(s - 3 W T) 

2 , . » 4 : ^ s * © *=6 2 

shifted 


0 . 7 ) 
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-aw T 

Hence if we define y - e 


^shifted “ . 


(5.8) 


The result of this transforaD.tion is th-?.t the poles/zeros are 
rotated by w T radians from its previous location. 

C 

If we now have the filter system function using the delay 
operator as 


H(z“^) = 


+ a,z“^ + 

o 1 

+ D,Z“^ + 

O 1 


-t gp 

H. z-5 


(5.9) 


the recursive relationship to implement this being 


( 3 . 10 ) 


where [x^] and [yj_] are the filter inputs and outputs respecti- 
vely, the conversion to a complex filter involves using new 
coefficients c^^g and such that [7J 


c^g = (5.11) 

dl (5.12) 

-This can be implemented recursively on a digital computer. 
Therefore, the design of a complex digital filter involves 

a) Design of a lowpass digital filter 

b) Shifting/ the filter characteristics to bandpass 
pbay n. cte-p-is'tics at the desired fre<|uea<^ . 
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3.5 Detc,ils of filters Designed 

A "banlc oi IP filters were used in the sinulc.tion. The 
IF frequency was 70 MHz cjid the RF bandwidth was 12 MHz. ¥e 
designed for Frequency Diversity Type I four third order 
Butterworth bandpass filters each having a bandwidth of 5 Mz 
as per details given in Sections 3.5 and 3.4. The lower two 
frequencies were assigned to mark and the upper two for space. 

The highest frequency to be passed being 76 ffiz, we had a sam- 
pling frequency of 152 MHz. 

For the study of frequency diversity Type II and Type III 
since we required a ba.nk of six IF filters in the assigned 
bandwidth of 12 MHz, we had to design new filters each with a 
passband of 2 MHz, centred at 65 Miz, 67 MHz, 69 I-Mz, 71 MHss, 73 
73 MHz, and 75 MHz. These were also third order Butterworth 
filters. The filters having centre frequencies of 65 I'lHz and 
75 Mz were assigned as diversity frequencies. The remaining 
four were assigned so that the lower two frequencies were meant 
for mark and the upper two frequencies were for space. 

3.6 Pseudo Han do m Hoise Generator 

In this simulation, complex Gaussian noise had to be 
added to the signal. The generation of complex Gaussian noise 
was done by generating pseudo random nxunbers which are uni- 
formly distributed and then using a transforma.tion to obtain 
a Gaussian distribution. 
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5 . 7 Greneration of Uniform Distr il autiori. 

The libro.ry subroutine for generation of pseudo randoia 
numbers uses the multiplicative congruentinl method. It 
consists of computing i +1"^^ random number using the relation- 
ship. 


X a (modulo m) (3.13) 

where x^ = i^^ pseudo random number, and 
a = a constant multiplier. 

The libraiy subroutine has a fixed starting seed. 

Therefore, the sequence of random numbers generated were the 
same for each simulation rxm since each simulation started 
with the same seed. Since the time for a single run was 
limited it was not possible to generate a continuous long 
sequence of random numbers. 

Since there was no provision for seeding the library 
subprogram for each run, it could not be used in this simu- 
lation. It was therefore essential to write a seperate function 
subprogram which could return the seed for subsequent iterations 
at the end of each run. Using these seeds as the starting 
seeds for the subsequent run, it was possible to generate 
suitable lengths of signal sequences without losing the 
generality of the entire sequence. 
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3 . 8 Generation of Graussian Distribution 

The Gaussian distribution was generated using the Box - 
Miller transformation 

s = (~2 logg r^)^ 0os(27ir2) (5.14) 

where r^^ and r 2 are two uniform random numbers between 0 and 1 
and s is the desired norna.1 distribution. 

To generate Gaussian distribution with a specified mean .■ 
and variance, we use the relationship 

X = ors + [1 (3.15) 

whore x is the sample from a normal distribution with specified 
mean u and variance a. s is as per eqn. (3.14). 

Since in this simulation we have to generate a complex 
Gaussian number, we generate two independent and identically 
distributed Gaussian numbers X and Y by eqns. (3.14) and (3.15) 
and form the complex sample 2 as given below 


2 = X + jY 


(3,16) 
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CHAPTER 4 

SYSTEM DESORIPTIOH AM DIVERSITY 

4.1 Sigml Plow Model 

The energy detection ESK system is given in [3] hns been 
used for this study. Pig. 4.1 shovrs the simplified block 
schematic of such a system. 

In the actual system the modixlated carrier through the 
channel tis passed throiigh a demoduli.tor . The output of the 
demodulator is fed through mark and space IP filters and then 
a decision is taken based upon which is the larger output. 

4 . 2 Simulated System 

The signal flow di-’gram used in the simulated system is 
shown in Pig. 4.2, 

In the case of the simulated system, we first form the 
channel using the tapped delay line model [1], [5]. The sign 
sequence generator generates the signal deciding a nark or a 
space based upon a uniform distribution. The frequency to be 
transmitted is chosen based upon firstly whether it is a mark 
or space and secondly the time slot in which it is to be 
transmitted. The signalling schemes that were used for 
various diversities are described later in this chapter. 

4.3 CSiannel Poxmilat i on 

The total number of taps , required in the delay line is 
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calcul-'ted as follows ; 


M TAPS = 


I 


-? . 

max '^nxn 


l/¥, 


1 laax min'^ 




where 


1 


A = 


max 


I'foltipafh delay 
Maximtim path delay 


= Mip.imum path delay 




= Sampling frequency 


(4.1) 


The delay power spectrum Q(6) was computed by evaluating 
the integral in eqn. (2.16) mentioned in Chapter 2 for various 
values of 6. The delay power spectrum is then sampled MTAPS 
times. These MTAPS volues of Q(6) are previously calculated 
and stored in the computer prior to the simulation. The 
complex tap gains are found by usiiig these MTAPS samples 
of the delay power spectrum Q(6) in the relationship 


^n ^^n 


Q(0,nA/) i DMSIvEPL 


1 < n ^ MTAPS 


where and b^ are Gaussian random variables of zero mean 
and unity varianee, UrSMPIi = the downssuapling factor, and 

¥ the samplir® frequen<^. 
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4.4 Modulation 

In tile modulator 1)0.361 upon the frequency, the signal 
is convolved with the channel using a recursive digital 
filtering technique described in [3]. ¥o use the relationship 

jwT jmwT 

^n=yn-*‘® ^n-1 - ® 

where the [y^] are sampled unit pulses, then the resulting 
[Xj^J will be equal to 

/ jwuT , 0 ^ n M-1 

XiCnt) = J (4.4) 

jo , n ^ M 

Thus we generate a M sample burst of a complex exponential 
with an angular frequency w. If now an arbitrary sampled 
sequence [yj_] is used as the input to eqn. (4.3) the resulting 
output will be the convolution of [y^J with the filter. 

PSE signals are rectangular bursts of complex exponentials. 
Therefore, if we use the complex tap gains of the channel as the 
set [y^] in eqn. (4.3) where w is the angular frequency of the 
PSE burst, then the output of the channel is the set [x^] 
as shown in eqn. (4.4). This relationship eqn. (4.5) is the 
shortest means for finding the effects of the channel on the 
signal and can be, for long channels, orders of magnitude 
faster than the usually used time domain convolution { 3 }, 
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4.5 Diiscital Filter 

The output of the modulator is now fed to digital If 
filters for the mark and space frev^uencies. These are third 
order Butterworth filters. These filters are complex, 
digital filters having a passband of width V a.bout some desired 
frequency _A.. . Unlike real filters, complex filters do not 
have passband centered around . 

These filters were designed using techniques given in 
[6], [7]. Ihe details of the digital filters design has 
air ready been given in Section 5.5. 

4.5 Gaussian Hoise. Generation 

Complex Gaussian Uoise of zero mean and unity variance 
is generated as described earlier in Section 3.8. This noise 
is filtered throijgh the respective IP filters and then is ' 
added to the signal after giving it necessary weightage to 
maintain a particular SITE. 

4 . 7 Energy Calculatio n 

The energy in the signal and noise, and the noise only 
are computed separately by summing the energy in the signal 
plus noise sequences and the noise sequences. This is done 
by summing the square of the magnitude of each sample. 

4.8 Decision 

A decision is arrived at this stage whether an error has 
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occured or not, This is done on the basis of comparing the 
energy contained in ©signal plus noise uitli the energy 
contained in only noise. An error has occured when the 
energy in noise filter is greater than the energy in the 
message carrying filter. Only under such circumstances, do 
we refer to the diversity channels, i.e. In case of a 
correct decision in the first channel we do not consider 
the energy of the diversity channel. 

4 . 9 Computation of Stat i stics 

On completion of a every 100 samples the number of errors 
that would have occured is printed out. This data is logged 
for each SNR so that the total number of errors can be 
computed for each SNR. Along with this the number of 

times the diversity channel was considered is also printed 
out . 

The 8-im of the simulation was to study the implicit 
diversity in the received waveform by virtue of frefjuency 
selective fading imposed by the channel on the received 
signal . 

4.11 Niversity Combining 

(Sie optiml method of combining the diversity is 

predetection iiaxiiaaih*atio oombining C4]» which involves 

l.l. }. h ^PUR 

CENIKAl USRAWf 
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■weighting each diversity channel received signal, hy a 
complex gain which is proportional to the conjugate of the 
channel gain and inversely proportional to the noise. However, 
in the case of this study due to paucity of computer time, we 
have used a modified form of equal gain post detection combining. 
In this case, we have referred to the diversity channel only 
in the case of an error. That is, if the first channel gives 

an error we refer to the diversity channel and sum up the 

energy of the signal and noise of the diversity channel and the 
first channel before arriving at a decision. 

Thus , if jj. 

2 C 21 (s-(t) + n.(t))^ - 51 (Hi (t))^3 0 

k=l i=.0 ^ i=0 ^ 

(4.5) 

where 

s(t) = signal amplitude 
n(t) = noise amplitude 

k = order of diversity 

M = number of samples per gate time, 
an error occurs. 

In our case, to reduce computer time we take a decision for 
k = 1 and in case of an error go for higher orders of diversity, 

4.12 Signalling Scheme 

The total allowable KH* bandwidth is divided into two 
equal, adjacent, overlapping frequencies called the 
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and ’space’ bands. Bach I'-lark and Space division is 
further divided into B equal divisions each numbered from 
1 to H in each band. 

The PSE transmission utilized to tx*aiismit at time tj, 

1 ^ I ^ N, a gated sine va.ve at either frequency MCI) or 
S(I) based upon whether a mark or a space is required to be 
transmitted. In the case of diversity transmission, signal 
formats as shown later in this chapter were transmitted. The 
pulse width T for the system is a constant. Also, 

3? t ”"^1 (4»6) 

Also, if the multipath spread is to be limited to 6 seconds, 

N(ti_j_i - tj) ^ 6 (4.7) 

If B is the Tnaxi m uTTi allowable bandwidth and ¥ is the band- 
max 

width of the pulse then, 

4.15 Signal Pornat Used 

Ihree types of signalling forraats were used. Each one 
consisted of a frequency time pattern, fhese are dealt with 
in detail below. 

4»15*1 Frequency Pi^^sity gype I 

Figure 4.3 illustrates the frefuen qy diversity !]^e I. 
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The bandwidth allocation here bein^ 4 we send the 
lower two frequency bands for Iferk and the upper two frequency 
bands for space. Hence let us say a ?ferk occurs then we 
send fj^ in time slot t^^ and the diversity channel is fg in 
time slot ± 2 * Similarly, fox* Space we send f^ in time slot 
t^ and f^ in time slot t 2 . Zlhe data rate in this case is 
^ b/secs with a bandwidth occupancy of 4 If we assume 

w^ as ¥, the data rate packing is bits/sec/Hz. 

In detection, time gates of duration = 2T are 
applied after the appropriate bandpass filters. Oomputation 
of the energy in the gated waveform provides the decision 
variables. Thus for the tJark filter time gate output woxild 
contain signal plus noise 

Wjjj(t) = S(t) + n(t) 

While the Space filter would contain only noise 
¥g(t) = n(t). 

Thus an error would occur if 

M-l< / 2 M>1 

51 - H 

i=0' i=0 

In the case of Type I, since we have a kind of second order 
diversity, the decision is taken after considering both the 
channels, an error occurs if 
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S^(t) + n^(t) 



0 


4 . 13,2 Frequency Diversity Type II 

Figure 4.4 illustrates the signalling format used in the 
diversity type II. 

In this case, we have two diversity channels at either 
end of the allocated bandwidth. The rest of the bandwidth 
is divided into four bands, the lower two frequency bands like 
earlier are for the mark and the upper two for the space frequ- 
encies. The diversity channel is chosen as the basis of whether 
a mark or space is sent and the time slot in which it is sent. 
For example, if a mark is sent in time slot t^, the diversity 
channel is So, the frequencies transmitted are f 2 and fg. 

Similarly, in time slot t 2 the frequencies are f^^ and f^ for 
mark. In the case of space, in time slot t^^ the frequencies 
Sent are f^ and fg and in time slot t 2 , the frequencies are 
fg and f-j^. 

The channel being the same for both the main channel and 
the diversity channel, the noise is constant during a time 
slot . 

Thus, for a mark in time slot tj^, we sum up the energy 
in mark filter fg and diversity filter fg arid compass?® it 
with the sum of the two noise energies in the respa^^re 
space filters. 
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Like earlier, an error occurs if 


2 M-1 

k=l i=0 


!i(t) 


+ n^(t) 


M-1 

r 

i=0 


n^(t) 


] 0 (4.11) 


The data rate in this case is l/T bits per second with a band- 
with allocation of 6 Again, assuming as ¥ the data 

rate packing factor in this case is bits/sec/Hsj. 

4.13.3 Frequency Diversity Type - III 

In this case the signal format used is as shown in Fig. 

4.5. 

This is similar to the frequency diversity case with the 
essential difference being that two time slots are used for 
each signal. Here, if ssy a iferk occurs, we transmit 
frequencies ±2 and fg in time slot t2_, f^ and f^^ in time 
slot ± 2 ' ' case of space, the frequencies are f^ and fg 

in time slot t^ and f^^ and fg in time slot tg. 

An error occurs in this case if 


4 M-1 

zc r 

k=l i=0 


S 


i(t) 


+ 


n^(t) 



(4.12) 


The data rate in this case is 1/2T bits per second, with 
a bandwidth allocation of 6 W^. Again assuming as ¥, 
the data rate packing factor in this case is 1/12T¥ bits/sec/Ez. 
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Pig. 4.1 Sigoal Plow Diagram 
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MARK INfRCV KROM OTHIR 



R6.4.2 SIGNAl FLOW DIAGRAM OF SIMULATED SYSTEM 
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CHAPTER 5 

SIMULATION AND RESULTS 


5.1 Introduction 

All the parameters for the simulation have been tahen from 
proposed or existing mobile troposcatter systems. A Fortran 
program for the event to event simulation of the ^stem was 
written. The program is modular in construction with the main 
program controlling all the subroutines. The system was 
simulated on the Digital Computer IBM 7044. 

5.2 Main Program 

The flow chart of the main program is given in Fig. 5.1. 
The details are as under : 

Block ’A* ; The seeds for the generation of the uniform random 

numbers for deciding mark/ space and generating a Gaussian 
distribution is read in as data. The weights for realizing 
the various signal-to-noise ratios is also fed in here. The 
centre frequencies for the various IF filters is also an input 
parameter. The total number of samples for each computer run 
NSAB5P and the SNRs for which the simulation is done is also 
decid ed. 

Block ; In this block the various counters/flags that are 

required for computing various statistics at the end of a 
simulation are initialized. These are for counting : 



(a) Humber of signal samples i.e, mark/space 

(b) Humber of errors 

The total energy in signal plus noise and the total 
energy in noise alone are also initialized to zero here, so 
as to evaluate the SHE. The various flags that are required 
during the simulation run are also initialized here. 

Block * G* : In this block the various storage locations that 
require initializing for each iteration are initialized. 

These include : 

(a) Signal plus noise energy of each signal sample 

(b) Hoise energy of each signal sample 

(c) Flag to decide signal sequence and frequency to he 
transmitted which is dependent upon type of diversity 
considered. 

(d) The iteration counter is incremented here 

Block 'D* : In the block D1 it is decided whether the signal 

is a mark or a space based upon a uniform distribution. In 
the second part of this block D2 on the basis of the si^al 
being a mark or a 8pace,the time slot and the type of diversity being 
used, the signal pattern that has to be sent through the system is 
is formulated. The system patterns are given in Pigs. 4.3, 4.4, 4.5 
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Block *E* : The channel is formed by calling the subroutine 

CHAHITEL. This consists of generating the 25 complex tap gains 
of the tapped delay. 

Block * P* ; Digital filtering of the complex tap gains is 
done prior to convolving the channel in subroutine DlfflPl, 

The Digital filtering operation is done at the particular IP 
filter centre frequency which has been decided earlier. The 
bandwidth of .the IF filter is 3 MHz in the Type I and 2 MHz 
in Type II and II J. 

Block * G* j The signal is convolved with the filtered channel 
using subroutine GONVOL here using the algorithm 

mentioned in Section 4.4 described earlier in Chapter IV. 

¥e use 26 samples of complex exponential at the specified 
frequency which are convolved with the 25 complex tap gains of 
the tapped delay line. The output of the convolution is 
50 samples of convolved frequency. 

Block : Complex Hoise is generated in this block using a 
subroutine HOIS that generates a complex Gaussian distributed 
number. These are then digitally filtered using subroutine 
DIGIFL. 

Block : This block reduces the noise as per the SHE 


required for the simulation using a weighting factor. 
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Block *K* ; In subroutine EiffiSGrT the suitably weighted fifty 

complex samples are individtially added to the fifty complex 
noise samples. The energy of signal plus noise is computed. 
Similarly the energy of noise is computed. This is done by 
evaluating the square of the magnitude in each case and summing 
it over 50 samples in each gate time. 

Block * L' i It is decided at this stage whether the total 
energy of the fifty signal plus noise samples is greater than 
total energy of fifty noise samples alone. 

Block : If in the first instance the signal plus noise 

energy is less than noise energy, the diversity channel is 
considered by branching to points 20,50 ami 40, which are 
predecided depending on the signalling format for which the 
system is being tested. The main channel is first tested 
and if the signal plus noise energy is greater than noise 
energy alone we take it as a no error case and skip to the 
next sample. 

An error occurs only vrhen the total signal plus 
noise energy is less than the total noise energy for all 
diversity channels. This is a variation over the eqxaal 
ratio post detection combining that should have been done in 
an actual system. This was done so as to reduce the ocm'- 
putation time. 
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Block * N* : On completion of the specified number of samples 

irSAMP» for the computer run, statistics are computed in this 
block. These statistics are : 

(a) Total number of iterations, i.e, Number of signal 
samples for which the system was tested NSAMP 

(b) Number of errors 

(c) Number of times diversity channel was used 

(d) Signal to noise ratio 

5.3 Results 

!I3ie system was simulated for signal -to-noise ratios 
(SNp) of 0,3,6,10,20,23 and 24.5 dB for types I and II, In 
Type III the signals-to-noise ratios for which the system was 
tested are as above except for 23 & 24.5 dB^he small signal-to- 
noise ratios of 0,5 an! 6 dB were simulated with a runs of 
over 1000 samples for each SNR. The larger SBRs of 10,20,23 
24.3 dB were tested with a simulation length of over 4000 

®h© graph showing probability of error versus signal~to- 
noise ratio are shown in Fig. 5.2. The probability of error 
is better in all types of diversity than the 2 PSK case as 
simulated in [8]. 

It is founi that Type I has a smaller probability of 
aS compared to Type II. This should be expected sinoa tfea 
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bandwidth of the IF filters is just 2 MHz in Type II as compared 
to 3 MHz in Type I, 

Frequency diversity Type III being a kind of quadruple 
diversity, it has the lowest probability of error, A compari- 
son of data rates shows that the diversity Type II has a better data 
rate than I and III for the same bandwidth. Among the three 
cases it is seen that probability of error is least in Type III 
even though there is a reduction in data rate by half as 
compared to Type II using the same bandwidth. 



START 


Read starting Seeds, SiiR weights 
IP centre freq.uencies, number ofj 
iterations, SliEs to be used. / 


BIOCE A 


Initialize counters for samples,' 
errors, flags, storage for total 
signal plus noise energy, Roise 
energy, SHE. 


bxoce: B 


Increment iteration counter, 
initialize signal plus noise 
energy, noise energy of each 
iteration. 


BLOCS C 


Generate Mark or Space by a uni- 
form distribution 


BLOCK D1 


Choose IP freq. to be used based 
on mark/space, diversity used and 
time slot 


BLOCK D2 
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Pig. 5.1 Plow Chart 
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CEiPTER 6 

coimusioi 

In this thesis we have studied some techniques of infcnnd 
diversity that would be suitable for a mobile troposcatter 
system. Incorporating such a technique of diversity in a 
system has the advantage of reduction of bulk that is 
essential for mobile troposcatter systems. The conventio n al 
forms of diversity like space and frequency do not permit such 
a drastic reduction in bulk though they do have better data 
rates. Inband frequenqy diversity also reduces alignment 
time as compared to conventional forms of diversity. 

Simulation of a traposcatter system like in this case takes 
many hours of computation time on a fa.st digital computer. 

In this case at various stn-ges we had to introduce modifications 
so as to reduce the ccaaputer time. With this in view, the 
type of combining had to be modified, so as to take a decision 
bathed ©a the main channel itself if correct and adding diversity 
channel ec^j^gjiee only in the case of errors. This will give 
quite an optimistic estimation of the probability of error, 
^■'urther, the simulation length had to be curtailed drastically 
due to prolonged computation time. It is necessary to have 
atleast eight to ten such simulation runs for each SNR before 
we can fix any probability of error with some degree of confi- 
dence. 
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The results obtained here give only an idea about the 
performance of the systems simulated but longer and more 
simulation lains are required to have a more accurate estimate 
of the probability of error, V/e expect better performance at 
higher SHRs while using inband frequency diversity. There is 
a drastic reduction in the data rate by these forms of diversity 
as compared to the conventional forms like frequency and space, 

A fourth order diversity scheme using a time frequency pattern 
like Type I involves reduction in the data rate to a quarter 
of the original data rate. However, with usage of codes as 
in !I^ype III, it should be possible to increase this data rate 
so that we can still get a fourth order diversity with a 
reduction in data rate by only 0.5. Many such codes can be 
used for transmission of signals. 

One such possibility is the usage of the signalling 
format shown in Pig. 6.1, This is another case of coded 
diversity. Here when a mark occurs, we transmit in time 
Blot frequencies f^ and f^ and in time slot tg frequencies 
fg aM whereas in the case of space we transmit frequencies 
fg and 'in time slot t^ and frequencies f^ and f^ in time 
slot tg. Since this requii'es only fovir bands, this has a 
higher data packing factor as compared to Type III mentioned 
earlier. It should also have a probability of error compara- 
ble to Type III. 



49 


One of the promising new areas of signal design for 
fading dispersive channels is the usage of error correction 
codes. D. Chase in [10] gives a combined coding and mnd ula - 
tion approach for communication over such channels. It is 
intuitively possible to visualize that considerable perfor- 
mance improvements is possible with coding without any loss 
in data rate beyond the normal redundancy of the code^ 
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